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Preface 


This  thesis  is  a  preliminary  design  study  of  a  toroidal  electro¬ 
static  analyzer  with  multiple,  off-center  slits.  The  use  of  off-center 
slits  provides  one  with  multiple  data  channels,  with  each  channel 
sampling  at  a  different  energy.  Thus,  energy  spectra  can  be  obtained 
in  a  fraction  of  the  time  normally  required  by  an  analyzer  with  a 
single  entrance-exit  slit  pair. 

I  would  like  to  thank  my  thesis  advisor.  Major  James  Lange,  for 
his  encouragement  and  frequent  help  which  was  so  essential  for  the 
completion  of  my  thesis.  1  am  also  grateful  to  lLt  John  Glessner  and 
lLt  Glenn  Param,  whose  conversations  were  of  frequent  aid. 

Finally,  I  wish  to  give  my  special  thanks  to  my  wife,  Carol,  for 
her  patience  and  assistance  throughout  my  thesis.  If  not  for  her 
loving  support,  this  paper  might  never  have  reached  its  completion. 


Robert  LaVerne  Hartley,  Jr. 
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Abstract 


The  design  of  a  toroidal,  electrostatic  analyzer  is  discussed 
with  special  emphasis  on  the  instrument  providing  multiple  channels 
of  data  by  sampling  different  points  of  the  energy  spectra 
simultaneously.  Two  additional  results  of  this  study  are  that  a 
toroidal  analyzer  can  have  four  times  the  resolution  of  a  spherical 
analyzer  and  can  cover  sixteen  times  as  much  of  the  energy  spectra 
(with  a  fixed  electrode  voltage). 

The  electrostatics,  equations  of  motion,  and  the  optic  equations 
are  discussed.  The  applications  chapter  discusses  the  use  of  multiple, 
off-center  slits,  the  minimization  of  the  size  of  the  instrument,  the 
trade-offs  that  must  be  made  in  designing  an  instrument,  the  effect  of 
field  errors,  the  minimization  of  the  effects  of  the  field  errors, 
and  the  total  beam  current  as  seen  by  the  detector.  Included  in  the 
Appendix  are  design  curves  that  should  be  of  assistance  to  the  reader 
in  designing  his  own  instrument. 


I.  Introduction 
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energy  ratio.  By  varying  the  voltage  applied  to  the  electrodes,  a 
complete  q/To  spectra  can  be  obtained.  Furthermore,  if  q  is  known, 
then  the  energy  spectra  can  be  measured  directly.  In  this  paper  the 
charge  will  be  folded  into  the  centripetal  force  and  will  not  show 
explicitly  as  q/T<j.  The  energy  will  be  given  by  itself  without 
reference  to  the  charge. 

The  four  most  common  types  of  electrostatic  analyzers  have 
parallel  plate  electrodes,  concentric  cylindrical  electrodes, 
concentric  spherical  electrodes,  and  toroidal  electrodes.  The 
toroidal  electrode  surfaces  are  described  by  two  distinct  radii, 
r  and  R,  with  different  centers  (see  Figure  2).  The  toroidal  elec¬ 
trodes  are  more  general  in  that  they  contain  the  cylindrical 
electrodes,  R  •»  «>,  and  the  spherical  electrodes,  R  =  r,  as  special 
cases . 

The  parallel  plate  analyzer  is  restricted  in  its  resolution  due  to 
the  limited  angle  through  which  it  can  deflect  the  beam.  The 
cylindrical,  spherical,  and  toroidal  analyzers  all  have  better 
resolution  than  the  parallel  plate  analyzer.  Because  of  their 
rotational  symmetry,  they  can  deflect  the  ion  beam  through  a  larger 
angle,  thus  achieving  better  energy  dispersion.  However,  the 
cylindrical  analyzer  has  no  axial  focusing,  thus  the  beam  spreads  out 
axially,  which  lowers  the  beam  intensity.  Due  to  their  double  curva¬ 
ture,  both  the  spherical  and  toroidal  analyzers  have  axial  focusing, 
which  helps  to  maintain  the  beam  intensity.  Since  the  cylindrical  and 
spherical  analyzers  are  both  special  cases  of  the  toroidal  analyzer, 
only  the  toroidal  analyzer  will  be  developed  in  this  paper.  When 


appropriate,  cylindrical  and  spherical  analyzer  formulas  will  be 
provided  by  taking  the  proper  limits  of  the  toroidal  analyzer  formulas. 

Electrostatic  analyzers  generally  have  high  resolution  and  high 
transmission  curves,  but  they  are  limited  to  an  energy  range  of  about 
10  ev  to  100  kev.  The  lower  energy  limit  occurs  when  the  electrode 
fields  are  reduced  to  about  the  same  magnitude  as  the  self-fields  of 
the  ion  beam.  The  upper  energy  limit  is  due  to  voltage  breakdown 
across  the  electrodes. 

Electrostatic  analyzers  are  much  lighter  and  have  simpler  power 
requirements  than  do  analyzers  which  use  magnets.  They  are  thus 
well-suited  for  space  applications,  such  as  solar  wind  experiments. 
They  are  also  useful  in  plasma  experiments  in  the  laboratory. 

B.  Problem 

Generally,  the  electrostatic  analyzer  provides  only  a  single  data 
channel.  In  certain  applications,  as  in  pulsed  plasma  experiments, 
there  is  only  a  very  limited,  total  observation  time.  This  makes  it 
difficult  to  gather  sufficient  data  in  a  single  experiment.  It  is 
hoped  that  by  extending  the  cross  radius  R  through  a  full  2it 
revolution,  that  is  the  radius  that  sweeps  perpendicular  to  the  main 
ion  path  (see  Figure  2),  that  an  instrument  can  be  designed  that 
provides  multiple  data  channels,  each  channel  sampling  a  different 
energy. 

C.  Coordinate  System 

The  coordinate  system  that  will  be  employed  is  the  cylindrical 

coordinate  system  r,  <J>,  z  shown  in  Figures  2  and  3.  The  radial  radii 

of  the  electrodes  are  r  and  r. .  The  cross  or  axial  radii  of  the 

a  b 


electrodes  are  R  and  R,  .  It  will  be  assumed  that  R  and  R.  have  a 
a  d  a  d 

common  center  of  curvature.  The  angular  extent  of  the  electrodes  in 
the  Z  *  constant  plane  is  $  ■  $e.  The  angular  extent  of  the  electrodes 
in  the  $  *  constant  plane  is  6,  where  6  will  be  allowed  to  sweep 
through  a  full  2n. 

The  mid-electrode  surface,  halfway  between  the  two  electrodes,  is 
described  by  the  radii  r  *  ag  and  R  *  Rg,  as  shown  in  Figure  3.  The 
central  path  through  the  electrodes  is  given  by  r  =  ag  and  z  =  0.  It 
will  be  assumed  that  the  ions  enter  the  electrodes  near  this  central 
path  with  small  entry  angles  such  that  the  ions  always  remain  in  the 
close  vicinity  of  the  central  path. 

In  most  places  in  this  paper,  the  dimensionless  coordinates  p  and 
£  will  be  used  where  p  and  £  are  defined  by 


(r  -  ae) 


(1-4) 


e 


(1-5) 


The  assumptions  in  the  previous  paragraph  require  that  p  and  £,  be 
ouch  less  than  one  at  all  times. 


D.  Approach 

The  electrostatic  equations  will  be  obtained  from  the  literature 
and  modified  for  use  in  this  paper.  Equations  of  motion  will  be  set 
up  and  solved,  thus  describing  the  trajectory  of  an  ion  through  the 
electrodes.  The  solutions  to  the  equations  of  motion  will  then  be 
manipulated  so  as  to  produce  optic  equations,  which  provide  a  more 
convenient  description  of  the  ion  behavior.  The  resolution  and 


transmission  will  also  be  developed  at  that  time.  Applications  of  the 
optic  equations  will  then  be  discussed  with  special  attention  paid  to 
the  behavior  of  off-center  slits.  It  is  the  use  of  off-center  slits 
that  permits  the  analyzer  to  sample  the  energy  spectra  at  various 
points  simultaneously,  thus  providing  multiple  channels  of  data. 

The  trade-offs  that  are  required  in  designing  an  actual  instrument 
will  then  be  discussed,  with  an  accent  on  minimizing  the  size  of  the 
instrument.  Design  curves  involving  the  various  system  parameters  are 
presented  in  the  Appendix  to  assist  the  reader  in  designing  an 
instrument  suitable  to  his  particular  needs. 

The  effects  of  errors  in  the  series  solutions  to  the  electric 
fields  are  discussed  with  respect  to  how  they  impact  on  the  resolution 
and  transmission.  The  effects  of  these  errors  can  be  minimized  with 
the  use  of  secondary  entrance  slits  which  restrict  the  entrance  angles. 

Finally,  a  qualitative  discussion  is  given  on  the  statistical 
problem  of  ion  flux  through  the  various  slits  and  the  electrodes  as  it 
relates  to  the  total  ion  current  seen  by  the  detector. 


II.  Electrostatics 


A.  Overview 

The  equations  for  the  radial  and  axial  electric  fields,  the 
electrostatic  potential,  and  the  equipotential  surfaces  are  listed, 
as  derived  by  N.  Svartholm  (Ref  1:196-197).  The  variables  and 
constants  involved  are  defined  and  derived.  The  connection  between 
the  voltage,  ±Uo,  on  the  electrodes  and  the  central  path  electric 
field  strength,  Eo,  is  derived.  Finally,  approximations  are  made  to 
arrive  at  the  final  form  of  these  equations  that  will  be  used  in 
subsequent  chapters. 

B.  Electrostatic  Equations 

As  mentioned  in  the  Introduction,  the  electric  field  strength 
equations  were  obtained  from  the  literature.  We  will  assume  that  the 
electric  field  and  the  electric  potential  are  both  zero  outside  of  the 
electrodes.  Fringing  fields  will  be  neglected  in  this  treatment.  It 
will  further  be  assumed  that  the  electrodes  are  charged  to  ±U0  and 
that  the  zero  equipotential  surface  is  midway  between  the  two 
electrodes,  for  small  electrode  gaps.  The  electric  field  is  assumed 
to  possess  cylindrical  symmetry  with  respect  to  the  z-axis  and  mirror 
symmetry  with  respect  to  the  z  ■  0  plane.  If  the  electric  field  is 
given  in  the  z  -  0  plane  near  r  »  ag  by  the  series  expansion 

E(p)  ■  E0[l  ♦  B iP  ♦  02p2  ♦  6jP3  ♦  .  .  .]  (2-1) 


where 


(2-2) 


then  according  to  N.  Svartholm  (Ref  1:196),  the  electric  potential  and 
the  electric  field  strength  in  an  annulus  surrounding  r  «=  ag,  is 
given  by 


U(p,;)  =  -E0ae[p  ♦  jBiP2  +  ^p3  ♦  JBsP1* 

-  1(1  +  Bi);2  +  j(l  -  Bi  -  2B2 )PC2 


-  j(l  -  Bi  ♦  B2  ♦  3B3)P2C2 


♦  ^(1  -  Bi  ♦  4B2  ♦  6B3)C“] 


Er(p.^)  -  E0[l  +  BiP  +  B2p2  +  BsP3 


(2-3) 


+  1(1  -  Bi  ♦  Ba  ♦  36s )P2C  +  |(1  -  Bi  ♦  4B2  +  6B3)C3] 


E  (p,£>  -  E0[-(l  +  Bi)C  +  (1  -  Bi  -  2B2)PC 


-  (1  -  Bi  +  B2  +3B3 )p2C 
+  1(1  -  Bj  -  4B2  ♦  6B3);3] 


(2-4) 


(2-3) 


where 


z 

*  “  — 
e 


(2-6) 


The  potential  has  been  set  up  for  U  ■  0,  when  r  ■  ag,  z  ■  0. 
When  the  potential  is  set  to 


U(p,<;)  »  -E0aeO 


(2-7) 
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where  a  is  a  constant,  then  the  equipotential  surface  with  this 
particular  potential  value  is  given  by  (Ref  1:197)  as 

p  -  o  -  5B102  ♦  (jBi 2  -  582)°*  +  |[1  ♦  81  -  (1  ♦  Bi2  -  2B2)o];2 

(2-8) 


The  constants  Bi  and  B2  are  given  by  H.  Ewald  and  H.  Liebl 
(Ref  2:873)  as 


8l  *  “(I  + 

Ke 

(2-9) 

82 

«  1  +  f*  +  a  2(1  +  R  * )/2R  2 

Ke  e  e  e 

(2-10) 

where 

R  '  2  ~| 
e  dr  I 

|r*ae 

|  z«0 

(2-11) 

In  our  case,  we 

have 

R  »  R+  (r  -  a  ) 
e  e 

(2-12) 

therefore. 

R  *  -  1 
e 

(2-13) 

This  gives  us 

•«  ■ 1  *  %  *  & 

(2-14) 

It  appears  that  Ewald* s  values  for  Bi  and  B2  were  arrived  at  by 
one  of  his  associates,  R.  Albrecht  (Ref:3),  who  used  a  general,  but 
difficult,  conformal  mapping  approach  involving  figures  of  revolution. 

A  simpler  approach  was  suggested  by  H.  Wollnik  (Ref  4:165-167),  for 
the  special  case  when  R^  and  R^  are  concentric,  as  is  the  case  here. 

By  symmetry,  one  can  assume  that  the  equipotential  surfaces  have  the 
same  shape  as  the  electrodes.  Conservation  of  electric  flux  can  be 
used  to  relate  the  electric  field,  E0,  for  r  ■  ag,  to  the  electric 


tr/. 


field,  E,  for  r  t  a^.  This  gives  us,  see  figure  2, 


where 


E( rd(j>)  (Rd0)  -  E0(a  d<{>)  (R  d0) 

e  e 


r  =  a  ( 1  +  p) 
e 


R  =  R  ♦  (r  -  a  ) 
e  e 


Therefore , 


E  -  Ec 


•MM 

_1 _  _ 1_ 

+  P>J[^  +  j 


Equation  (2-18)  can  be  written,  to  second  order  in  p,  as 


E  *  E0  1(1  -  p  ♦  p2) 


(l  -  &>  ♦  S£>*) 


E  -  E0  1  ♦ 


This  gives  us  immediately 


•('  •  *)' 

•  « 


_  i  &e  2 

P  +  1  +  5=  ♦  p^ 

Ke  “e‘ 


-('  •  *) 

(l  *  S  *  tf) 


(2-15) 


(2-16) 


(2-17) 


(2-18) 


(2-19) 


(2-20) 


(2-21) 


(2-9) 


(2-14) 


These  values  for  and  $2  are  consistent  with  those  given  by 
Ewald.  It  should  be  noted  equations  (2-3),  (2-4),  and  (2-5)  apply  to 
either  a  cylindrical,  spherical,  or  toroidal  system.  It  is  the  choice 
of  ae  and  Re  that  determines  to  which  system  the  equations  apply. 

For  the  cylindrical  system  Re  ■  for  the  spherical  system  Re  =  a  , 


J 


first  order  in  p  or  £.  Applying  these  constraints  to  equations  (2-3), 
(2-4),  (2-5)  leaves  us  with 

U(p,0  =»  -E0ae  p  +  jBiP2  -  |(1  +  MS2  (2-26) 

Ef(p)  *  E o [ 1  +  Sip]  (2-27) 

E z(rJ  =  E0 [-( 1  +  Si);]  (2-28) 

It  is  these  three  equations  that  will  be  used  in  the  mechanics 
chapter  in  the  equations  of  motion. 


III.  Mechanics 


A.  Overview 

The  pattern  of  development  of  chapter  III  closely  follows  that 
used  by  H.  E.  Duckworth  and  S.  N.  Ghoshal  (Ref  5:210-221),  who  will 
not  be  further  cited.  Certain  notational  changes  are  made  to  improve 
the  similarity  in  form  between  the  radial  and  axial  equations  and  to 
generalize  the  radial  solution.  The  most  important  change  from 
Duckworth's  notation  is  the  transformation  of  the  radial  solution  to 
explicitly  include  the  ion  energy. 

The  primary  goal  of  the  mechanics  chapter  is  to  arrive  at  solutions 
for  the  vadial  and  axial  motion  of  the  ions  inside  of  the  electrodes. 
The  equations  of  motion  are  the  standard,  cylindrical  equations  of 
motion. 

By  invoking  conservation  of  energy,  the  centripetal  force  equation 
and  the  electrostatic  potential,  the  velocity  change  as  the  ion 
penetrates  the  field  boundary  is  derived.  This  result  is  coupled  with 
conservation  of  angular  momentum  to  get  a  solution  for  the  angular 
equation  of  motion. 

By  combining  the  solution  to  the  angular  equation  of  motion, 
centripetal  force  equation,  and  the  radial  electrostatic  field,  the 
radial  equation  of  motion  can  be  solved.  It  is  in  this  section  that  a 
connection  is  made  between  certain  math  constants  in  the  radial 
equation  and  the  energy  of  the  ions,  relative  to  the  mean  energy,  T0, 
to  which  the  system  is  tuned. 


15 


By  combining  the  solution  to  the  angular  equation,  the  centripetal 
force  equation,  and  the  axial  electric  field,  the  axial  equation  of 
motion  can  be  solved,  thus  completing  the  mechanics  chapter. 

B .  Equations  of  Motion 

The  equations  of  motion  are 

Mr  -  Mr<$>2  =  eEr  (3-1) 

M  £  (r2$)  =  eE^  =  0  (3-2) 

Mz  =  eEz  (3-3) 

-19 

The  change  of  the  ion  is  e  =  +1.609  x  10  coulomb,  and  the  ion  mass, 
M,  is  defined  to  be 


M  =  M0( 1+Y)  (3-4) 

where  Y  is  assumed  to  be  a  small  number.  M0  will  be  defined  shortly^by 
the  centripetal  force  equation. 

C .  Velocity  Change  at  Field  Boundary 

For  an  ion  of  mass,  M,  and  a  velocity 

v'  =  v0(l+B)  (3-5) 

outside  the  electric  field,  with  3  assumed  to  be  a  small  number,  the 
velocity,  v,  of  the  ion  just  after  it  enters  the  field  is  given  by  the 
energy  balance  equation, 

|  M( v ' 2 )  »  i  Mv?  +  eU  (3-6) 
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This  can  be  rearranged  as 


9 


’  * 


(3-7) 


For  an  ion  of  mass  Mo  and  velocity  vo  undergoing  pure  angular 
motion  in  the  field  along  p  *  o  and  z  =  0,  the  centripetal  force 
equation  is 


(3-8) 


This  gives 


eEoae  =  -MoVo2  (3-9) 

By  combining  equations  (2-9),  (2-26),  (3-4),  (3-5),  (3-7),  (3-9), 
rearranging,  and  expanding  to  second  order  in  p,  C,  6*  or  "f ,  we  arrive 
at 


This  is  the  velocity,  v,  of  an  ion  in  the  field  that  started  with  a 
velocity,  v'  *  Vo(l+8),  outside  of  the  field.  Equation  (3-10)  will  also 
be  needed  to  first  order: 


v  =  vo(l+B-p)  (3-11) 

D.  Angular  Equation  of  Motion 

The  ion  velocity  component  along  the  central  path  is  the  angular 
velocity,  v^,  which  is  given  by 

v^  «  vct.i  ;  (3-12) 
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where  0}  is  the  angle  between  the  ion  path  and  the  central  path,  v  is 

the  ion  velocity  inside  the  field  after  it  enters  at  the  point  (pi,£i). 

Since  the  entrance  angles  were  assumed  to  be  small,  we  can  use  the 

1  2 

approximation,  cos(Oi)  =  1  -  j  Oi  .  By  combining  this  with  the 
equations  for  v,  (3-10),  and  for  v^,  (3-12),  expanding  and  retaining 
only  second  order  terms  in  8.  Y»  Pi*  and  £i,  we  have 

v4>  -  vo£l+8-Pi+^  (Pi2-Ci2)+Pi(|+^-|  Oi2j  (3-13) 

We  can  also  write 


•  _  v  v 

<Pl  “  r  i  “  ae(  1+Pi ) 

We  thus  have 


(3-14) 


♦  i  *  JJ  ^l+8-2p i- j8p i+Yp i  +2p i 2+  5  (f^)  Cpi2-Ci  2)-|  o,2J  (3-15) 

The  second  equation  of  motion  gives  us  conservation  of  angular 
momentum,  from  which  we  can  write 


r2$  =  r  1 2<J  i  (3-16) 

By  substituting  the  equation  for  <j>i,  (3-15),  into  equation  (3-16)  and 
solving  for  <J>  to  second  order  in  8*  Y,  p,  pj,  £j,  and  oj,  we  arrive  at 

♦  -  ^l+8-2p-2p8+3pz-p j 2+j8P i+PiY-j  o,2+i  (|&)  (Pl2-C,2)j  (3-17) 

This  equation  describes  the  angular  motion  of  an  ion,  but  it  is 
nonlinear  and  coupled.  The  <£  equation  was  derived  to  second  order  only 
to  keep  from  losing  first  order  terms.  To  first  order,  the  ^equation  is 


$  -  ( i+8-2p) 

ae 


(3-18) 
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This  equation  for  $  will  be  used  in  the  radial  equation  of  motion. 


It  will  also  be  useful  to  drop  the  $  equation  to  zeroth  order  and 
integrate  with  respect  to  time. 

4>  -  ~  (3-19) 

ae 

This  equation  will  be  useful  later  in  eliminating  (t),  from  the  radial 
and  axial  solutions,  in  favor  of  ($). 

E.  Radial  Equation  of  Motion 

As  r  •  ae(l+p),  we  have  r  «  aep;  therefore,  the  radial  equation  of 
motion  (3-1)  can  be  transformed  to 


Ma.p  -  Mae(l+p)  (<|>)2  «  eEp 


(3-20) 


By  combining  the  equations  for  Bi  (2-9),  Ep  (2-27),  M  (3-4),  the 
centripetal  force  (3-9),  and  £  (3-18),  with  the  radial  equation  (3-20), 
we  arrive,  to  first  order  in  8,  Y,  and  p,  at 


B  .  (£)’  [<Y,26)  -  p(2-g)] 


By  defining  the  following  three  new  constants 


e2  -  *£ 
Re 


X2  *  2  -  e2 
r  Y  +  23 


(3-21) 


(3-22) 

(3-23) 

(3-24) 


we  can  transform  the  radial  equation  of  motion  to  its  final  form 
shown  below 


U 


p(<t>)  »  A2C08(X4>)  ♦  B2sin(X<t>)  ♦  6  (3-30) 

For  X2  <  0,  we  have 

p C <t> )  *  A3cosh(X<fc)  +  B3*sinh(X4>)  ♦  6  (3-31) 

Since  X2  < 0,  X  is  purely  imaginary.  We  can  then  define  X*  as  follows, 
with  i  »  J-=r . 

iX*  -  i  ,/c!-2'  -  X  0-32) 

By  defining  Bj  «  iB3*,  the  third  solution  becomes 

p($)  ■  A3Cos(X*0)  ♦  B3sin(X*4>)  ♦  6  (3-33) 
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By  redefining  X2  to  be 


X2  -  2  -  eJ 


(3-34) 


then  X*  -  X,  and  the  third  solution  becomes  identical  in  form  to  the 


second  solution. 


The  boundary  conditions  are  declared  to  be 


P<0)  -  Pi 


(3-35) 


(3-36) 


where  ar'  will  be  defined  later  in  the  optics  chapter  as  the  radial 
entry  angle.  Applying  these  conditions  gives  us 


X2  «  0  p(4>)  -  -®r'4>  ♦  Pi 


(3-37) 


X2  f  0  p(0)  *  (Pi-6)cos(X4>)  — sin(X<j>)  +  6  (3-28) 

In  the  limit  as  X2-*’0,  the  solution  for  X2  +  0  reduces  to  the 
X2  *  0  solution.  It  is  thus  seen  that  the  trig  solution  (3-38) 
contains  both  the  linear  solution  and  the  hyperbolic  solution. 

The  parameter  6  can  be  related  to  the  kinetic  energy  of  the  ion. 
We  can  define  T  to  be  the  energy  of  an  ion  outside  the  field  with 
mass  M  and  velocity  v'f  and  To  to  be  the  kinetic  energy  of  an  ion  of 
mass  Mo  and  velocity  vo .  Therefore,  we  have 


T  -  |m(v* ) 2  -  |mo(1+Y)vo2(1+S)Z 


T0(l+Y)  (1+B): 


(3-39) 
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To  first  order  in  &  and  Y,  we  have 


T  -  T0(l  ♦  26  ♦  Y) 

By  applying  the  definition  for  6  (3-24),  we  have 

T  -  T0(l  ♦  X26)  (3-40) 


Solving  for  6  gives  us 


(3-41) 


By  defining  the  relative  energy,  ri»  as 


(3-42) 


we  have 

(3-43) 

It  should  be  noted  that  T  is  the  unperturbed  kinetic  energy  of  an 
ion  outside  of  the  electrostatic  field,  whereas  T0  is  the  kinetic 
energy  of  an  ion  inside  the  field.  Specifically,  T0  is  the  kinetic 
energy  to  which  the  analyzer  is  tuned.  It  is  this  energy  that  will 
permit  an  ion  to  track  the  central  path  through  the  electrodes. 

By  substituting  this  expression  for  6  (3-42)  into  equation  (3-38), 
we  arrive  at  the  final  form  of  the  radial  solution: 

p  ■  (pi  -  cos(X4>)  -  sin(X4>)  ♦  ^  (3-44) 

This  form  of  the  radial  solution  is  in  terms  of  the  relative  ion 
energy,  the  point  of  entry,  the  entry  angle,  and  the  electrode 
geometry. 


6  * 
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F.  Axial  Equation  of  Motion 


As  z  ■  agC,  we  have  z  ■  agc;  therefore,  the  axial  equation  of 
motion  (3-3)  can  be  transformed  to 


Mae£  ■  eE^ 


(3-45) 


By  combining  the  equations  for  Bi  (2-9),  Ej.  (2-28),  M  (3-4),  and 
the  centripetal  force  (3-9),  with  the  axial  equation  (3-44),  we 
arrive,  to  first  order  in  8,  Y,  and  C,  at 


5  -  o 


(3-46) 


The  standard  solution  to  this  equation  is 


C(t)  ■  Acos^e^^-j  ♦  Bsin^e—^ 


By  invoking  the  equation  for  $  (3-19),  we  have 


(3-47) 


C(<J>)  «  Acos(e4>)  ♦  Bsin(e<(>) 


(3-48) 


The  boundary  conditions  are  declared  to  be 


C(0)  -  Cl 


(3-49) 


(3-50) 


where  az'  will  be  defined  later  in  the  optics  chapter  as  the  axial 
entry  angle.  Applying  these  conditions,  we  arrive  at 


C(<t>)  *  Cic°s(c4>)  -  -&-sin(e<|>) 


(3-51) 


This  form  of  the  axial  solution  is  in  terms  of  the  point  of  entry, 
the  entry  angle,  and  the  electrode  geometry.  It  differs  in  form  from 
the  radial  solution  only  in  that  the  relative  energy  does  not  appear 
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in  the  axial  solution.  This  difference  is  principally  connected  to 
the  fact  that  the  axial  electric  field  is  much  smaller  than  the 
radial  electric  field. 

G.  Summary  of  Key  Mechanics  Equations 

For  the  convenience  of  the  reader,  certain  key  mechanics 
equations  which  will  be  needed  in  the  optics  chapter  are  listed  here 
for  future  reference. 


e2  .  is. 

Re 

(3-22) 

X2  -  |2  -  e2 | 

(3-34) 

_  T  -  To 
n  To 

(3-42) 

t? 

p(<J>)  *  ^pi  -  cos(X4>)  sin(X4>)  ♦  ^jr 

(3-44) 

C(<t>)  -  Cicos(e<|))  -  sin(e4>) 

(3-51) 

A.  Overview 


The  pattern  of  development  of  chapter  IV  closely  follows  that  used 
by  H.  E.  Duckworth  and  S.  N.  Ghoshal  (Ref  5:210-221),  who  will  not  be 
further  cited.  Many  of  the  toroidal  equations  presented  in  this 
chapter  were  developed  in  analogy  to  the  cylindrical  equations 
produced  by  Duckworth  and  Ghoshal.  Certain  notational  changes, 
consistent  with  chapter  III,  are  made  to  improve  the  similarity  in 
form  between  the  radial  and  axial  equations  and  to  permit  the  equations 
to  explicitly  contain  the  relative  ion  energy. 

In  this  chapter,  coordinate  systems  are  established  in  the  regions 
outside  the  entrance  and  exit  to  the  electrodes.  The  behavior  of  the 
ions  in  these  regions  is  worked  out  and  refraction  at  the  field 
boundary  is  considered. 

By  combining  the  behavior  of  the  ions  outside  the  electrodes  with 
the  solutions  to  the  radial  and  axial  equations  of  motion  inside  the 
electrodes,  the  locations  of  the  radial  and  axial  image  planes  are 
determined. 

Further  manipulation  of  the  equations  which  specify  the  locations 
of  the  image  planes  yields  the  location  of  the  radial  and  axial  focal 
points,  the  radial  and  axial  "focal  lengths,"  and  the  radial  and  axial 
Newtonian  lens  equations.  It  must  be  remembered  that,  while  these 
equations  are  patterned  after  the  thin  optical  lens  equations,  the 
electrostatic  analyzer  is  not  a  thin  optical  lens  and  thus, some  of  the 
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dx'  ■  -aed<J>.  The  positive  y'-direction  matches  the  positive  r  or 
p-direction;  thus,  at  the  field  boundary  dy'  =  dr  =>  agdp.  The  positive 
z'-direction  matches  the  positive  z  or  £ -direction;  thus,  at  the  field 
boundary  dz*  »  dz  ■  aed£. 

The  angle  ar  is  defined  to  be  the  angle  between  the  x'-axis  and 
the  projection  of  the  ion  path  onto  the  x'y'-plane  (see  Figure  4).  It 
should  be  noted  that  if  >  0,  evaluated  at  the  field  boundary,  then 
p  is  decreasing.  This  is  because  the  ion  is  moving  in  the  (-x') 
direction. 

The  angle  «z  is  defined  to  be  the  angle  between  the  x'-axis  and  the 
projection  of  the  ion  path  onto  the  x'z'-plane  (see  Figure  4).  It 
should  be  noted  that  if  <xz  >  0,  evaluated  at  the  field  boundary,  then 
4  is  decreasing.  As  in  the  radial  case,  this  is  because  the  ion  is 
moving  in  the  (-x')  direction. 

The  angles  and  will  both  be  assumed  to  be  small  angles. 

This  is  required  in  order  to  be  consistent  with  the  earlier  assumption 
in  the  electrostatics  chapter  that  the  ion  path  must  remain  in  the 
close  vicinity  of  the  central  path  through  the  electrodes.  With  the 
assumption  of  small  angles,  we  can  write 
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Figure  5  -  Entrance  Angles 


By  the  use  of  simple  geometry,  we  can  arrive  at  expressions  for 
and  cc^  in  terms  of  the  coordinates.  Let  us  assume  that  the  ion 
originated  from  the  off-axis  point  A1  with  the  coordinates  x'  =2.', 
y'  ■  aePo»  z'  =>  ae£o  (see  Figure  5).  Let  us  further  assume  that  the 
ion  enters  the  field  at  a  point  with  the  coordinates  x'  *  0,  y'  *  aepj 
z'  =*  ae^i .  From  the  geometry,  it  is  apparent  that 

tan  (“^)  =  (p0  -  Pj)  (4-3) 

tan  (°^  -  «^  «  (Co  “  Cl)  (*“*) 

Equations  (4-3)  and  (4-4)  will  permit  the  eventual  elimination  of 

<*  and  a  in  favor  of  the  coordinates  of  the  ion  source  and  the 
r  z 

coordinates  of  the  ion's  point  of  entry  into  the  electric  field. 

C.  Exit  Coordinate  System 

The  exit  coordinate  system  x",  y",  z"  has  its  origin  located  at 
P  ■  0,  c  ■  0*  4>  *  4>e  (see  Figure  4).  The  positive  x"-direction  points 
away  from  the  electrodes,  but  this  time  it  is  in  the  direction  of  the 
ion  motion;  thus,  at  the  origin  dx"  = 


+a  d4>.  The  positive  y"  and 


z"-directions  match  up  with  the  positive  p  and  c -directions 
respectively;  thus,  at  the  field  boundary  dy"  *  aedp  an(*  dz"  =  aed£. 

The  angle  is  defined  to  be  the  angle  between  the  x"-axis  and 
the  projection  of  the  ion  path  onto  the  x"y"-plane.  As  the  ion  is  now 

II 

moving  in  the  direction  of  increasing  x" ,  the  angle  «  >  0  corresponds 

to  increasing  p. 

The  angle  «"  is  defined  to  be  the  angle  between  the  x"-axis  and  the 
z  • 

projection  of  the  ion  path  onto  the  x"z"-plane.  We  now  have  (similar 
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to  the  radial  case)  >  0  corresponding  to  increasing 

We  will  assume  that  and  are  small  angles.  This  is  consistent 

i  i 

with  the  assumption  that  and  cc^  are  small  angles  and  with  time 
reversal  symmetry.  With  the  assumption  of  small  angles,  we  can  write 


dp  .  <&;; 

d<j>  dx '  „  . 

x"*0 


M  II 

tan  (oc  )  =  cc 
r  r 


(4-5) 


dr  I  dz" 


x"=0 


»  ii 

tan  (“  )  -  cc 
z  z 


(4-6) 


By  using  the  equations  for  p  (3-44)  and  for  c  (3-51),  we  can 
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evaluate  the  dif ferentials  in  expressions  for  a  and  a  as 


cc"  =  *£ 

r  d$ 


-X^Pi  -  sin(X<t>e)  -  ar  cos(X4>e)  (4-7) 


d$ 


-cc i  sin(c<t>e)  -  “  cos(e<{>  ) 


(4-8) 


<t>=4>€ 
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This  gives  us  and  in  terms  of  the  entry  angles,  the  electrode 


geometry,  and  the  coordinates  of  the  ion's  point  of  entry  into  the 
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electric  field.  We  will  need  and  «  shortly  to  compute  the  ion 
motion  after  it  exits  from  the  electrodes. 

D.  Refraction  at  Field  Boundary 

It  would  be  apropos  to  investigate  the  effect  of  ion  path  bending 
or  refraction  at  the  field  boundary  at  this  point  before  continuing. 


For  an  ion  penetrating  a  boundary  with  a  finite  potential  change, 
the  velocity  component  parallel  to  the  boundary,  VM ,  is  essentially 
unchanged.  It  is  the  velocity  component  perpendicular  to  the  boundary, 
V^,  that  is  affected  by  the  change  in  potential  energy  (Ref  5:173). 

We  thus  have  (see  Figure  6)  outside  the  boundary,  sin  (a^)  =  V||/V^Q, 
and  inside  the  boundary,  sin  (<^)  »  V(|/V^^.  Solving  these  two  equations 
for  V 1 1  and  equating  them  gives  us  an  electrostatic  Snell's  law: 

V^0  sin  («Q)  =  V  £  sin  («^)  (4-9) 
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or  for  small  angles 


V  a  =  V  oc . 
IP  O  ±1  1 


(4-10) 


We  can  equate  Vj^  to  V'  (3-5),  a0  to  <*r  (4-3),  and  to  V  (3-11) 
with  p  »  Pj.  We  therefore  have 


V'oc'  =  V*.  (4-11) 

r  l 

or 

V0(l  ♦  8)  (PG  -  Pi)  =  V0(l  ♦  6  -  Pi)«i  (4-12) 

Solving  for  oc^  to  first  order  in  8,  PQ  and  Pi,we  have 

ai  *  fr  <Po  ~  Pi) 

0 

»  «*  (4-13) 

i. 

I 

By  exact  analogy,  we  also  have  in  the  axial  case,  •  oc^,  where 
<x^  is  now  an  axial  angle  instead  of  a  radial  angle.  Thus  to  first 
order,  we  may  ignore  any  ion  refraction  at  the  field  boundary.  This 
is  true  at  both  the  entrance  or  the  exit. 

Refraction  could  not  be  ignored  in  any  study  which  was  interested 
in  examining  higher-order  effects,  such  as  image  aberrations.  To  get 
some  feel  for  the  effect  of  the  higher  order  terras,  we  start  again 
with  the  electrostatic  Snell’s  law  (4-9).  By  assuming  that  the  angles 
are  small,  we  can  replace  the  sine  function  in  equation  (4-9)  with  the 
tangent  function: 


0^' 


w<* 


Of  greater  interest  is  the  effect  of  refraction  on  the  image 
location.  By  applying  an  equation  that  will  be  developed  later,  the 
radial  Newtonian  lens  equation  (gf  ■  0,  2”  *  2.'),  we  get 


2r"t  .  _ (1  ♦  3) _ 

r  ^1  +  3  -  Pi  +  Pi2  +  |  PiB) 


(4-16) 


i'J  is  the  radial  image  location  that  will  be  derived  later, 
assuming  no  refraction.  2r"  is  the  radial  image  location  when  refrac¬ 
tion  is  considered.  Figure  9,  on  the  previous  page,  shows  equation  (4-16) 
graphed  for  several  values  of  3*  The  ion  energy  T  is  greater  than/less 
than  T0  as  3  is  greater  than/less  than  zero. 

Figure  9  shows  that  ions  with  an  energy  T  <  To  tend  to  focus  in  a 
shorter  distance  than  those  ions  with  T  >  To.  Figure  9  also  shows  that 
ions  with  an  entry  point  such  that  p!  >  0  are  focused  farther  out  than 
those  ions  with  Pi  <0.  Just  exactly  where  an  ion  ends  up  depends  on 
its  particular  values  of  T  and  pi.  In  general,  as  the  ion  energy 
deviates  from  To  and  the  entry  point  varies  from  pi  »  0;  the  ion  image 
changes  to  a  new  location.  If  the  exit  slit  is  located  at  2^',  then  the 
refraction  effects  lead  to  a  loss  of  ions.  Along  Pi  =  0,  the  error  in 
the  image  location  is  slightly  stronger  for  T  <  To  than  it  is  for 
T  >  To.  These  effects  distort  both  the  transmission  curve  and  the 
resolution.  While  these  effects  are  definitely  of  interest,  for  purposes 
of  this  study,  they  will  be  ignored  in  the  rest  of  this  paper. 


or  for  small  angles 


n/»\  ti  •* 

y  (x  )  *  y  (0)  +  x  ar 


(4-18) 


At  x  =0,  we  can  equate  y  (0)  with  aep(<{>e).  We  can  thus  combine 
the  above  equation  (4-18)  with  the  equation  for  p  (3-44)  evaluated  at 
4>  *  4>e»  equation  for  «”  (4-7),  and  the  equation  for  (4-3),  to 


y"(x")  »  n  g  cos(X$e)  +  ~  sin(X4>e) 


-P°  sin(X4»e)  +  cos(X<}>e) 


♦Pi  ae(l  +  jcos(X<J>e)  +  (^JT  ~  x"xjsin(X<t>e)  (4-19) 

It  is  desirable  that  all  of  the  ions  of  a  given  energy,  that 
originate  from  one  given  point  A*,  be  focused  at  the  same  point  a|! 
after  they  exit  the  electrodes.  This  is  possible  if  and  only  if  the 
coefficient  of  pi  in  equation  (4-19)  vanishes.  By  setting  the 
coefficient  of  pi  to  zero,  we  can  solve  it  for  x",  thus  determining 
the  location  l"  of  the  radial  image  plane  in  terms  of  the  other  system 
parameters.  This  gives  us 


i,  _  B "  ae2  aeft1  .  v 

x"  =  2-r  -  -CT-  ♦  v”  cot(x<t>e) 


S.'  -  4^  cot(X<M 


(4-20) 


F.  Axial  Image  Location 


Similar  to  the  radial  case,  the  Z"  coordinate  of  an  ion  after  it 


exits  the  electrodes  is  given  by 


Z"(X")  ■  Z"(0)  +  X"  tan(«") 


(4-21) 
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or  for  small  angles 


Z"(X")  -  Z"(0)  ♦  X"  (4-22) 

At  X"  -  0,  we  can  equate  Z"(0)  with  ae£(<t>e)«  We  can  thus  combine 
the  above  equation  (4-22)  with  the  equation  for  £  (3-51)  evaluated 
at  $  ■  $e»  the  equation  for  «z  (4-8),  and  the  equation  for  (4-4), 
to  get 


As  in  the  radial  case,  we  desire  that  all  of  the  ions  of  a  given 
energy,  .that  originate  from  one  given  point  A',  be  focused  at  the  same 
point  Az  after  they  exit  the  electrodes.  This  is  possible  if  and  only 
if  the  coefficient  of  in  equation  (4-23)  vanishes.  By  setting  the 
coefficient  of  to  zero,  we  can  solve  it  for  x",  thus  determining  the 
location  SL |J  of  the  axial  image  plane  in  terms  of  the  other  system 
parameters.  This  gives  us 


V  -  cot(e$e)  (4-24) 

Note  the  extreme  similarity  of  form  of  the  axial  equation  (4-24)  with 
the  radial  equation  (4-20).  This  will  occur  often  as  we  proceed  with 
the  Optics  chapter. 


G.  Radial  and  Axial  Focal  Poincs 


c 


It  is  of  interest  to  know  where  a  parallel  beam  of  ions  into  the 
system  will  be  focused.  This  is  easily  obtained  by  taking  the  limit 
of  and  £ z  as  £'  goes  to  infinity.  From  equations  (4—20)  and  (4-24), 


we  have 


g”  ~  £im  £p  =  4^-  cot ( X4>e ) 

n  I  _ _  ^ 


(4-25) 


gz  =  £im  £z  *  cot(e<t>e) 
£'-h» 


(4-26) 


II  IV 

The  locations  of  the  radial  and  the  axial  focal  points  are  gf  and  gz 
respectively. 

It  is  also  of  interest  to  know  what  £'  would  have  to  be  in  order  to 
provide  for  a  parallel  beam  of  ions  out  of  the  system.  In  this  case 

II  II  ..... 

£r  and  £z  separately  would  be  equal  to  infinity.  Solving  equations 
(4-20)  and  (4-24)  for  £'  (now  called  £r  and  £z  respectively),  we  get 


£r  COt(X<$>e) 


£r  -  cot(X<f>e ) 


(4-27) 


£z  =  *§T  +  cot(etpe) 


£z  -  cot(e4>e) 


(4-28) 


Taking  the  limits  as  £r  and  £z'  go  to  infinity  gives  us 


g'  =  £im  £*  ■  cot(x<J>e) 

1  0  ••  L  A 

-XX) 


(4-29) 


g'  =  £im  £z  -  ^  cot(eif>e) 
"s  n  n  .  e  c 

z 


(4-30) 
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As  expected,  for  reasons  of  symmetry,  we  have  gr  =  gr  =  gr  and 
gz  *  g2  5  gz.  Thus  the  focal  points  (measured  from  the  ends  of  the 
electrodes)  are  located  at 


gr  ■  ^  cot(X$e) 


gz  *  “  cot(e<t>e) 


(4-31 

(4-32 


H .  Newtonian  Lens  Equations 

Substituting  gr  (4-31)  and  gz  (4-32)  into  the  equations  for 
£r'  (4-20)  and  (4-24)  gives  us 

2 


r  -  *  *'Sr 


(4-33 


‘  Sr 


a 


(4-34 


**  -  S2 


The  above  equations  can  be  rearranged  as 


ar  -  gr)  a*  -  gr)  -  Zfr  +  gr: 


a2  -  8z)  (£•  -  gz>  -  +  g; 


(4-35 

(4-36 


These  two  equations  are  almost  identical  in  form  with  a  Newtonian 
lens  equation.  If  the  right-hand  sides  are  defined  to  be  ff2  and  fz2 
respectively,  then  the  form  would  match.  Thus,  the  Newtonian  lens 
equations  are  defined  to  be 


with 


2  2 
r  2  -  ae  .  2  _ *e _ 

fr  -  "fr  gr  Xzsin*(X<t>e ) 

,  2  -  2  2  _ ae 2 _ 

z  ~  e2  +  ®z  “  eisin^c(pe) 


(4-39) 

(4-40) 


It  must  be  noted  that  while  fr  and  fz  are  called  the  radial  and 
axial  focal  lengths,  they  are  merely  definitions  of  convenience  that 
permit  us  to  write  the  lens  equations,  which  are  themselves  only  a 
convenient  form.  It  is  gr  and  gz  that  are  actually  important. 


I .  Post-Electrode  Ion  Location 

The  location  of  an  ion  is  given  by  y"(x")  (4-19)  and  z"(x")  (4-23). 
It  will  be  convenient  later  if  these  two  equations  are  transformed 
now  to  simpler  forms. 

II 

By  setting  x"  «  £r,  the  px  coefficient  vanishes  in  equation  (4-19), 
giving  us 


y"Ur")  -  ae 


$ 


-Po 


-cos(X$e)  +  1  +  X2.r  sin(X4>e) 


ae  sin(X4>e)  +  cos(X4'e) 

yt’  t* 


(4-41) 


The  radial,  Newtonian  lens  equation  (4-37)  can  be  rewritten  as 


T^r 


Sr 


+  Sr 


(4-42) 


By  combining  the  equation  for  y"(2.^  (4-41)  with  the  equations  for 
2-r  (4-42),  gr  (4-31),  and  fr  (4-39),  we  can  obtain  a  new  expression  for 


n 

[l  +  /  \1 

n  /  I  fr  \ 

Y7 

• 

(4-43) 
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This  equation  expresses  the  radial  location  o£  an  ion  in  the  radial 
image  plane  as  a  function  of  the  ion's  relative  energy,  n,  and  its 
radial  starting  position,  po. 

By  setting  x"  *  fi.” ,  the  £1  coefficient  vanishes  in  equation  (4-23), 


giving  us 


3  r  sin(e<J>e)  ♦  l'z  cos(c<{>0) 

■  ~  -ae4o  e£i  ^T 


(4-44) 


The  axial,  Newtonian  lens  equation  (4-38)  can  be  rewritten  as 


‘'z  =  £1  _  „  +  8z 


(4-45) 


By  combining  the  equation  for  Z"(iz)  (4-44)  with  the  equations  for 
lz  (4-45),  gz  (4-32),  and  fz  (4-40),  we  can  obtain  a  new  expression 
for  z"(C): 


Z  ( iz )  -  -agCo^,  *  gJ 


(4-46) 


This  equation  expresses  the  axial  location  of  an  ion  in  the  axial 
image  plane  as  a  function  of  the  ion's  axial  starting  position,  C,o 
Note  that  it  is  independent  of  the  ion's  relative  energy,  n* 


J .  Radial  and  Axial  Magnification 

Let  the  entrance  slit  have  a  radial  width  of  wr  *  2aep0, 
symmetrically  located  about  p  =  0,  i.e.,  the  edges  of  the  slit  are 
located  at  p  ■  ±p0.  The  width  of  the  image  of  the  entrance  slit  in  the 
radial  image  plane  is 


wr  -  y"<Po)  -  y"(-Po> 


(4-47) 
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By  equation  (4-43),  we  have 


*r  =  ae  * 


-ae  7T 


(4-48) 


■  -U'°‘(l^Tr)  ‘  (l^) 


(4-49) 


Wr 

If  we  define  the  radial  magnification  to  be  Gr  =  «-4*  ,  we  then  have 


=r  -  (r^hr) 


(4-50) 


Similarly,  let  the  entrance  slit  have  an  axial  height  of 
w^  =  2aeCoi  symmetrically  located  about  £  =  0,  i.e.,  the  edges  of  the 
slit  are  located  at  £  *  ±^0 .  The  height  of  the  image  of  the  entrance 
slit  in  the  axial  image  plane  is 


w"  -  z"Uo)  -  Z"(-Co) 


(4-51) 


By  equation  (4-46),  we  have 


“z  -  (f  f--g-2) 

"*  ■  -fe)  ’  -«  (f^i;) 


(4-52) 


(4-53) 


If  we  define  the  axial  magnification  to  be  Gz ^  — f-  ,  we  then  have 


(4-54) 


It  is  now  possible  to  write  the  post-electrode  ion  location 
equations  (4-43)  and  (4-46)  in  their  most  compact  form.  Substituting 


the  equations  (4-50)  and  (4-54)  for  Gr  and  Gz,  respectively,  into 
equations  (4-43)  and  (4-46)  yields 


.  n  [1  +  Gr]  -  p0G 

y  (n,Po)  -  ae  jf  r  K  r 


(4-55) 


Z"(Co)  =  -aeCoG, 


(4-56) 


It  must  be  noted  that,  in  this  section,  Po  and  £q  are  performing 
double  duty.  In  the  derivation  of  Gr  and  Gz,  po  and  £0  were  specific, 
fixed  values.  Tn  the  two  equations  above,  they  are  again  arbitrary, 
as  per  their  usual  usage. 


K.  Energy  Dispersion 

The  energy  is  a  measure  of  how  far  the  image  of  a  beam  of  ions  is 
displaced  because  the  ion  energy  is  not  the  energy  to  which  the  system 
is  currently  tuned.  This  displacement  is  called  D. 

The  radial  energy  dispersion,  Dr,  and  the  radial  energy  dispersion 
coefficient,  Kr,  are  given  by  (Ref  4:179),  (Ref  5:215),  (Ref  6:27) 


Dr  =  Krn 


(4-57) 


where 


.  ay” 


(4-58) 


With  the  aid  of  the  equation  for  y"  (4-55),  we  have 


Kr  “  x?  (1  +  Gr) 


(4-59) 


Similarly,  the  axial  energy  dispersion,  Dz ,  and  the  axial 
dispersion  coefficient,  Kz,  are  given  by 


m 


where 


(4-61) 


rr 


However,  with  the  help  of  the  equation  for  z"  (4-56),  we  have 

Kz  =  0  (4-62) 

Therefore,  Dz  =  0  (4-63) 

There  is  no  axial  energy  dispersion,  at  least  not  to  first  order. 
This  is  understandable  as  the  major  velocity  components  lie  in  the 
radial  plane  and  the  axial  velocity  component  is  small.  One  could  try 
to  rework  the  axial  equations  as  a  second  order  approximation,  but  the 
axial  equation  of  motion  is  of  a  form  for  which  no  general  method  of 
solution  exists.  Therefore,  it  will  not  be  pursued  here. 

i 

It  should  be  noted  that  n  contains  To  in  the  denominator  and  thus, 
Dr  varies  inversely  with  the  mean  energy,  T0 . 

L .  Resolution 

The  absolute  resolution  is  the  reciprocal  of  the  minimum  difference 
in  absolute  energy  that  can  be  resolved  as  two  separate  lines  by  the 
system  (Ref  4:179).  The  relative  resolution  is  the  reciprocal  of  the 
minimum  difference  in  the  relative  energy  that  can  be  resolved  as  two 
separate  lines  when  the  system  is  tuned  for  the  absolute  energy,  T0. 

The  beam  image  moves  inward  or  outward  linearly  with  n,  as  h  is 
decreased  or  increased.  This  is  expressed  by  the  radial  dispersion 
equation 

Dr  =  K  n  =  yf-  (1  +  Gr)n  (4-57) 

X  (4-59) 
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By  inspection  of  the  y"  equation  (4-55) 


"(n,Po>  =  ag  {  7-7  (1  +  Gr)  -  PocJ 


(4-55) 


>  .  ■  ,  , 

it  can  be  seen  that,  for  a  positive  change  br  in  p0,  the  image  moves 


inward  by  br,  where 


br  -  -brGr 


(4-64) 


If  the  entrance  and  exit  slits  are  symmetrically  located  about 
p  =  0,  ( y *  =  y"  =  0),  and  have  widths  w'  and  w"  related  by 


w"  =  Grw' 


(4-65) 


then  a  beam  of  ions  of  energy  t)  =  0,  with  p0  =  0  and  “r  =  0,  will  have 

its  image  located  at  y"  *  0  (see  Figure  10).  If  the  relative  energy 

is  increased,  the  beam  image  will  move  outward  until  it  just  touches 

W**  GrW  * 

the  edge  of  the  exit  slit  located  at  y"  ■  —  *  — —  ♦  By  increasing 
w  * 

Po  to  - — ,  the  beam  image  will  move  back  toward  the  center  of  the  exit 
slit.  The  relative  energy  can  then  be  increased  again  until  it 
reaches  a  value  defined  as  p+  in  which  the  beam  just  touches  the  outer 
edge  of  the  entrance  slit  and  the  beam  image  just  touches  the  outer 
edge  of  the  exit  slit.  The  relative  energy  n+  is  the  maximum  relative 
energy  a  parallel  beam  of  ions  into  the  system  can  have  and  still  be 
transmitted  through  the  system. 

Similarly,  one  can  decrease  the  relative  energy  until  the  beam 

image  just  touches  the  inner  edge  of  the  exit  slit  located  at 

— W *  — GtW *  — w *  • 

y  *  — —  =  — — —  .  By  decreasing  po  to  -s —  ,  the  beam  image  will  move 

back  toward  the  center  of  the  exit  slit.  The  relative  energy  can  then 

be  decreased  further  until  it  reaches  a  value  defined  as  n  in  which 
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Figure  10  -  Graphic  Definition  of  n  and  n 
the  maximum/minimum  ion  energy 
transmittable  when  the  system  is 
tuned  to  r)  =  0,  T  *  To* 


the  beam  just  touches  the  inner  edge  of  the  entrance  slit  and  its 
image  just  touches  the  inner  edge  of  the  exit  slit.  The  relative 
energy  n  is  the  minimum  relative  energy  a  parallel  beam  of  ions  into 

the  system  can  have  and  still  be  transmitted  through  the  system. 

W**  GtV  *  w  *  + 

By  setting  y"  *»  -y-  =  — ^ —  ,  p0  =  yy  ,  and  n  =  n  ,  and  inserting 

these  values  into  the  equation  for  y"  (4-55),  we  get 


(4-66) 


Solving  this  equation  for  n  yields 


X2Grw* 
ae^Gr  +  1> 


(4-67) 


— w**  -Grw  *  — w  *  - 

By  setting  y"  =  — y-  =  — ^ —  ,  po  =  yy  ,  and  n  =  n  ,  and  inserting 

these  values  into  the  equation  for  y"  (4-55),  we  get 


-Grw'  .  rf  „  x  .  w'  „ 

2  -  ae  JT  <l  *  Gr>  ♦  257  r 


(4-68) 


Solving  this  equation  for  n  yields 


-X2Grw' 

ae(Gr  +  1) 


(4-69) 


The  relative  resolution  is  thus 


n+  n“  -  2X  Grw' 

n  ■ n  ■  Tjey  v  ii 


(4-70) 


(L  +  Gr ) 


(4-71) 


By  equation  (3-42) 


T  -  Tp 

To 


(3-42) 
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Therefore 


T+  =  T0n+  ♦  To 


T  a  Ton  +  To 

The  absolute  resolution  is  thus 


(4-72) 

(4-73) 


T+  -  T~  =  T0(n+  -  rf )  =  -~- 

Rrel 


(4-74) 


Therefore 


8,bS  ■  Vfo?  (l  -  £) 


r  ,  _  Rrel 

Rabs 


(4-75) 

(4-76) 


That  the  absolute  resolution  varies  inversely  with  T0  is 
understandable  in  that,  as  the  energy  of  the  ion  beam  goes  up,  it  is 
increasingly  difficult  to  separate  the  T+  ions  from  the  T  ions  by  the 
width  of  the  exit  slit.  This  result  agrees  favorably  with  Paolini's 
results  (Ref  7:584)  once  one  realizes  that  Paolini  has  written  the 
reciprocal  of  the  resolution. 

A  comparison  can  be  made  now  between  the  resolutions  of  the  three 
curved  systems:  cylindrical,  spherical,  and  toroidal.  For  reasons 
that  will  be  explained  in  the  Applications  chapter,  Gr  will  be  set 
equal  to  one.  By  taking  the  appropriate  values  of  Re,  we  get 


Cylindrical : 

Rrel  ‘  & 

(4-77) 

Spherical : 

Rrel  "  5* 

(4-78) 

Toroidal: 

Rrel  =  xftr 

(4-79) 
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The  spherical  case  has  twice  the  resolution  of  the  cylindrical 


case.  It  will  be  shown  in  the  Applications  chapter  that,  X2  ^  7. 


Therefore,  the  toroidal  case  has  as  much  as  four  times  the  resolution 
of  the  spherical  case.  This  is  clear  motivation  for  picking  the 
toroidal  system  over  the  other  systems. 


M.  Transmission 


The  transmission,  T(n)*is  the  fraction  of  a  beam,  at  a  specific 
energy  r],  that  passes  through  the  exit  slit  divided  by  the  fraction 
that  passed  through  the  entrance  slit.  We  can  ignore  ion  losses  inside 
the  electrodes  as  we  have  already  assumed  that  all  of  the  ions  remain 
near  the  central  path.  The  transmission  is  thus  the  fraction  overlap 
of  the  beam  image  with  the  exit  slit. 

To  get  the  transmission  we  need  to  know  the  location  of  the  edges 
of  the  beam  image.  The  edges  can  be  located  with  the  aid  of  the 
equation  for  y"  (4-55),  keeping  in  mind  that  there  is  an  image 


•  •  w 

inversion,  that  the  entrance  slit  edges  are  located  at  p0  *  ±  ■= — 

zae  t 


and  w"  ■  Grw'.  The  inner  edge  of  the  image,  y",  is  given  by 

la 


yL 


(1  *  Cr>  -  °n 


(4-80) 


The  outer  edge  of  the  image,  y^,  is  given  by 


yU  “  a*\w  (1  +  Gr)  +  IT  Gr 


(4-81) 


When  n  >  0,  yy  is  outside  the  area  of  the  exit  slit. 


Thus,  the 


transmission  is 

T(n) 


2  "  yL 


n  >  0 


(4-82) 
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or  T(n)  =  1  +  (1  +  ~)  n  n  >  o  (4-85) 

These  two  forms  can  be  combined  into  a  single  equation  that  is 
good  for  all  r|,  as  shown  below 

T(n)  =  1  "  (  l+~)|n|  (4-86) 

Recalling  the  equation  for  the  relative  resolution  q  (4-71),  we  can 
write  the  transmission  in  its  final  form  as 

T(q)  *1-2  Rrel  |q|  (4-87 ) 

A  comparison  is  readily  made  between  the  transmissions  of  the  three 
curved  systems.  As  developed  here,  all  three  have  transmission 
curves  that  have  the  shape  of  an  isosceles  triangle  with  a  peak  value 
of  100Z.  The  difference  between  the  curves  is  how  sharp  they  are,  and 
that  is  strictly  a  function  of  the  resolution.  Thus,  the  toroidal 
system  has  the  sharpest  transmission  curve. 


N .  Summary  of  Key  Optics  Equations 


rr 


For  Che  convenience  of  the  reader,  Che  key  optics  equations 
listed  here  for  future  reference. 

I 

Radial  Entrance  Angle,  «r 
“r  -  ff-  (Po  "  Pi) 

Axial  Entrance  Angle,  «z 

•  ae  /  v 

«Z  *  p-  (Co  -  Cl) 

Radial  Exit  Angle,  “j.' 

“r  *  -X^pi  -  -prj  sin(X<t>e)-“^  cos(X<t>e) 

,  .  II 

Axial  Exit  Angle,  “z 
<*z  «  -e  z, i  sin(e4>e)  -  a'z  cos(e<t>e) 

Location  of  Radial  Focal  Point,  gr 
gr  ■  cot(X4>e) 

Location  of  Axial  Focal  Point,  gz 
gz  =  —  cot  ( €4>e ) 

Radial  Newtonian  Lens  Equation 
<*r  -  8r>  <*'  -  gr)  -  f/ 

Axial  Newtonian  Lens  Equation 
(*z  "  gz)  <*'  -  8Z)  "  fz2 


are 


(4-3) 


(4-4) 


(4-7) 


(4-8) 


(4-31) 


(4-32) 


(4-37) 


(4-38) 
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Axial  Magnification,  Gz 


Radial  Image  Displacement,  b^"  (when  entrance  slit  moved  by  b^') 


b  "  «  -b  'G 
r  r  r 


(4-64) 


Radial,  Relative  Energy  Resolution,  Rre^ 
Rrel  =  2xfr  (X  +  G^) 


(4-71) 


Radial,  Absolute  Energy  Resolution,  Ra^s 


abs  Tc 


(4-76) 


Transmission,  T(n) 


T(n)  =  l  -  2Rrel  |n| 


(4-87) 
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V.  Applications 


A.  Overview 

The  results  from  the  previous  chapters  are  used  in  this  chapter  to 
work  out  the  behavior  of  off-center  slits.  It  is  the  use  of  off-center 
slits  that  gives  this  analyzer  its  chief  advantage  -  that  of  being  able 
to  sample  different  parts  of  the  energy  spectra  at  the  same  time  with 
a  single  voltage  on  the  electrodes. 

Criteria  are  developed  for  minimizing  tf  'i  size  of  the  instrument. 
This  leads  to  a  simplified  set  of  optic  equations  and  lays  the 
foundation  for  discussing  the  various  trade-offs  that  must  be  made  in 
designing  an  electrostatic  analyzer. 

Errors  in  the  series  solutions  are  then  considered  as  functions  of 
the  deviation  from  the  central  path.  A  discussion  is  given  on  how 
these  errors  affect  the  resolution  and  the  transmission.  This  leads  to 
choosing  secondary  entrance  slits  so  as  to  restrict  the  entry  angles 
so  that  the  ions  will  remain  in  regions  with  no  more  than  a  given 
error  in  the  series  solutions  to  the  electric  fields. 

Finally,  a  qualitative  discussion  is  given  on  how  to  determine  the 
total  beam  current  from  what  is  known  about  the  original  ion  population. 
This  information  is  needed  in  order  to  make  a  choice  of  the  entrance 
slit  size  when  trading-off  between  the  resolution  and  the  detector 
sensitivity.  It  is  also  needed  for  analyzing  the  data  provided  by  the 
instrument . 
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B.  Off-Center  Slits 


o 


One  of  the  attractions  of  a  toroidal,  electrostatic  analyzer  is 
that  the  radii,  and  R^,  can  be  allowed  to  sweep  through  a  full  2 tt 
(see  Figure  2).  This  eliminates  one  pair  of  edges  and  their  fringe 
fields.  Also,  ions  can  be  admitted  at  any  point  on  the  circle,  R  =  Rg , 
and  all  of  the  previous  results  still  hold,  although  a  new  coordinate 
system  is  needed  for  each  new  entry  point. 

These  new  coordinate  systems  differ  only  in  that  they  have  been 

rotated  about  the  center  of  curvature  for  R  (see  Figure  11).  This 

e 

rotation  yields  a  new  entry  point  and  a  new  coordinate  system,  but  the 
equations  used  for  the  new  coordinate  system  are  identical  in  form  to 
the  original  equations  used  for  the  original  entry  point  and  the 
original  coordinate  system. 

The  advantage  of  the  additional  entry  points  is  that  their  entrance 
and  exit  slits  can  be  placed  at  various  off-center  locations.  This 
permits  the  user  to  sample  different  mean  energies  simultaneously. 

As  the  voltage  on  the  electrodes  is  changed,  the  multiple  channels 
sweep  different  segments  of  the  energy  spectra  at  the  same  time. 

This  results  in  a  complete  spectra  in  far  less  time,  which  can  be 
critical  in  experiments  with  a  limited  total  observation  time.  A 
specific  example  that  comes  to  mind  is  a  pulsed  plasma  device. 

The  equation  which  describes  the  behavior  of  off-center  slits  is 

y"  «  ae{X2n(l  +  Gf )  -  PoGr)  (4-55) 

This  can  be  rewritten  as 
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By  letting  po  be  the  center  of  an  off-center  entrance  slit  and  y"  be 
the  center  of  the  off-center  exit  slit,  equation  (5-1)  determines  the 
mean,  relative  energy  that  is  being  sampled  by  that  pair  of  slits. 

The  voltage  that  is  being  applied  to  the  electrodes  is  hidden 
inside  of  n«  The  relative  energy,  n*  is  given  by 


(3-42) 


1  2 
where  T0  =  -  M0v 

From  the  centripetal  force  equation  (3-9),  we  have 


1 

2 


M0v02 


-eEoae 

2 


From  the  electrostatics  section,  we  have 


(3-9) 


E0 


where 


(2-24) 


(2-25) 


where  d  is  the  electrode  separation.  Therefore,  we  have 


(5-2) 


Equation  (5-2)  gives  the  mean,  absolute  energy.  To,  for  which  the 
on-center  slits  of  the  analyzer  are  tuned,  in  terms  of  the  electrode 
geometry,  the  unit  charge,  and  the  applied  voltage,  Uq.  It  must  be 
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remembered  that  Eo  is  intrinsically  negative;  therefore,  equation  (5-2) 
yields  a  positive  value  for  To. 


By  combining  the  equations  for  T  (5-1),  To  (5-2),  and  n  (3-42),  we 
arrive  at 


This  expression  gives  the  mean,  absolute  energy,  T,  for  an  off-center 
pair  of  slits,  in  terms  of  the  location  of  the  entrance  and  exit  slits, 
the  unit  charge,  the  electrode  geometry,  and  the  voltage  applied  to 
the  electrodes. 

The  range  of  energies  that  can  be  covered, for  a  fixed  voltage,  is 
determined  by  the  available  range  of  values  for  y"  and  p0.  To  avoid 
ambiguity,  for  each  radial,  R,  there  should  be  only  one  entrance  slit. 
The  number  of  exit  slits  is  limited  only  by  the  size  of  the  detector 
behind  each  slit.  The  pattern  of  exit  slits  can  be  repeated  over  and 
over  again,  each  group  with  a  different  entrance  slit  and  thus 
sampling  a  different  segment  of  the  energy  spectrum.  The  number  of 
groups  is  limited  only  by  the  size  of  the  detectors. 

The  absolute  energy  T  (5-3)  has  been  plotted  as  T/eUo  versus  e2 
in  the  Appendix  (see  Figures  13,  14,  15),  where  e  is  defined  to  be  the 


58 


unit  charge.  In  these  three  figures,  I  have  set  the  radial 
magnification  to  G  ■  1;  the  electrode  gap  to  d  =  ae/50;  and  the 

•  •  V** 

maximum  displacement  of  the  entrance/exit  slits  to  *—  =  Po  =  ±.01. 

ae 

In  Figure  13,  the  upper  energy  limit,  Tu/eUo,  of  the  off-center  slits 

v" 

is  plotted,  so  we  have  ^—  =  Po  =  +.01.  In  Figure  14,  the  lower  energy 

ae 

limit,  T^/eUo,  of  the  off-center  slits  is  plotted,  so  we  have 

y" 

=  po  =  -.01.  In  Figure  15,  the  difference  between  Figures  13  and 
ae 

14  has  been  plotted  as  AT/eUo.  This  difference  is  the  maximum  range 
of  energies  (divided  by  ello)  that  the  system  can  sample  simultaneously 
at  a  fixed  voltage.  For  a  given  value  of  e2,  the  energy,  To/eUo, 
being  sampled  by  the  on-center  slits  is  the  average  of  the  corre¬ 
sponding  values  from  Figures  13  and  14. 

To  use  these  figures:  first,  choose  a  value  for  e2  and  then  read 
T  Til  A  T 

— off  the  given  curves;  second,  multiply  by  a 
particular  value  U0  for  the  electrode  voltage.  This  will  give  you 
T  ,  T  ,  and  AT  in  terms  of  ev.  T0 ,  in  ev,  is  the  average  of  Tu  and  , 

To  illustrate  this  let  us  pick  e2  *  1.94.  We  then  have  (in  ev) , 

T^  =  58U0,  T^  =  42Uo,  and  AT  *  16U0.  For  Uo  *  160  volts,  we  then  get 
T^  =  9.28  kev,  To  =  8.0  kev,  T^  *  6.72  kev.  AT  is  2.56  kev,  or  about 
322  of  T0. 

By  comparison,  for  a  cylindrical  system,  we  have  ez  =  0.  Therefore, 
in  ev,  T^  =  50.2Uo,  T^  =  49.7U0,  and  AT  =  .5U0.  AT  is  about  12  of  T0. 

For  a  spherical  system,  e2  *  1.  Therefore,  in  ev,  we  have  T^  =  5O.5U0, 
T^  *  49.5U0,  and  AT  *  Uo.  AT  is  about  2%  of  T0.  For  a  toroidal  system, 
the  optimum  value  of  e2  is  1.94  (this  is  discussed  later  in  this  chapter) 


This  gives  us,  in  ev,  Tu  =  58Uo,  T^  =  42Uo,  and  AT  =  16U0.  AT  is 
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about  32Z  of  To.  The  toroidal  system  clearly  has  the  advantage  over 
the  other  two  systems. 

The  energy  range  sampled,  AT,  grows  as  |e2|  ■+  2.  By  good  fortune, 
the  resolution  also  is  maximized  by  ( e 2 [  •*  2  (to  be  shown  later).  It 
is  the  improvement  in  the  energy  range  and  resolution  that  makes  the 
toroidal  system  more  attractive  than  the  spherical  system. 

A  discrete  sampling  of  the  energy  spectra  can  be  obtained  by 
driving  the  applied  voltage  with  a  staircase  waveform.  If  the  steps 
have  a  slight  ramp  shape,  then  a  continuous  spectra  can  be  obtained 
while  the  jumps  prevent  unnecessary  overlap. 

Due  to  the  linearity  of  the  equations  involved,  the  resolution  of 
an  off-center  pair  of  slits  is  the  same  as  it  is  for  the  on-center 


slits,  and  is  given  by 


R  ■■ 

rel  2X2w 


*  ('^) 


(4-71) 


The  transmission  curve  is  unchanged  in  shape,  but  it  is  shifted 
in  the  relative  energy.  If  n1  is  the  relative  energy  for  which  the 
off-center  pair  of  slits  is  tuned,  as  per  equation  (5-1),  then  the 
transmission  is  given  by 


T(n)  -  l  -  2  Rrel ( I n |  -  n’) 


(5-4) 


C.  Instrument  Size-Maximum  Compactness 


If  one  neglects  the  size  of  the  instrument  for  a  moment  and  tries 
instead  to  maximize  the  resolution,  one  quickly  discovers  that  no 
physically  realizable  maximum  exists  with  respect  to  the  parameters 
ag  and  Rg.  A  maximum  does  exist  for  the  parameter  in  a  form  which 
requires  X  3  0;  therefore,  X$e  *  0*  which  forces  the  slits  outwards  to 
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infinity.  If  X,',  the  object  distance,  is  set  to  X,'  =  K*gr,  a  multiple 
of  the  distance  to  the  focal  point,  the  problem  remains  for  the 
parameters  ag  and  Rg,  while  the  maximum  with  respect  to  <f>e  now 
requires  X4>e  =  ntl,  n  =  0,2,4  ....  Unfortunately,  this  choice  of  X<t>e 
requires  that  the  entrance  and  exit  slits  be  located  at  infinity. 

The  above  discussion  immediately  brings  home  the  need  to  hold  down 
on  the  size  of  the  instrument,  especially  the  separation  of  the  two 
slits.  The  Newtonian  lens  equation  gives  us 

Ur"  -  gr)  (£'  -  gr)  -  fr2  (4-37) 

It  can  be  shown  with  elementary  calculus  that  if  the  product 

AB  =  constant,  then  the  sum,  A  +  B,  is  a  minimum  if  A  =  B.  This  gives 

us  the  requirement 

SLr"  *  SI'  (5-5) 

Therefore,  equation  (4-37)  becomes 

U'  -  g  )2  =  f  2  (5-6) 

r  r 

It  is  obvious  that  ('  is  a  minimum  if 

gr  =  0  (5-7) 

At  first  thought,  it  would  seem  that  having  the  focal  point 
located  at  zero  is  unreasonable,  but  it  must  be  remembered  that  this 
is  not  a  thin  system.  It  is  entirely  possible  to  have  the  focus 
located  at  zero.  Having  the  focal  point  located  at  zero  means  that  it 
is  located  in  the  plane  that  contains  the  exit  end  of  the  electrodes. 


By  equation  (4-31),  we  have 


nr 


g  «■  4r-  cot(X4>  )  (4-31) 

r  a  e 

gr  ■  0,  implies  that 

-  *  (5-8) 

e  2 

We  could  set.  equal  to  higher,  odd  multiples  of  ^  ,  but  we  must 

have  less  than  2tt  and  this  would  require  that  X  increase. 

Unfortunately,  increasing  X  decreases  the  resolution.  Thus,  the  best 

resolution  is  maintained  by  discarding  the  higher,  odd  multiples  of 

Therefore,  the  best  compactness  can  be  accomplished  by  requiring  that 

SL  "  =  £'  and  X4>  =  ^  . 

r  e  2 

D.  Simplified  Optic  Equations 

The  requirements  imposed  by  equations  (5-5)  and  (5-8)  (from  the 
previous  section)  give  us  the  following  simplified  optic  equations: 


*e 


(X2  «  |2  -  e2 | ) 

(‘2  "  ft) 


g  =  0 


o'  *  0  "  »  f  .  £fiL 

r  X 


-ill-*  g 

ae  6z 

T  "  8z 


f  = 


ae 


‘"(i) 


es 


(5-9) 

(5-10) 

(5-11) 

(5-12) 

(5-13) 

(5-14) 


T(n)  *  i  -j^fr  (111  -  n‘>  (5-18) 

(for  on-center  slits  n *  ■  0) 


The  last  two  equations  can  be  further  simplified  if  we  set 

—2 

w*  «=  Ka„.  Typically,  K  has  values  on  the  order  of  10  .  This  gives  us 


rel  “  W 


T(n)  ■  i  -  v4r  (|n|  -  n') 


(5-19) 


(5-20) 


These  simplified  equations  are  graphically  represented  in  the 
Appendix.  In  several  cases,  the  equation  has  been  divided  by  a£  before 
plotting  so  as  to  get  a  curve  of  greater  utility.  The  curves  have  been 
plotted  with  the  restriction  that  <pe  <  2tt  and  most  of  them  have  a 
common  e2  horizontal  axis.  These  curves  should  make  it  easier  for  the 
reader  to  select  parameters  for  his  own  instrument. 


E.  Trade-offs  and  Constraints 


The  selectable  parameters  are  ag ,  Rg ,  <J>  ,  K  (the  slit  width),  and 
d  (the  electrode  gap).  These  parameters  determine  £',  £r"  (the  slit 
locations),  Rfe^  (the  resolution),  and  strongly  affect  the  voltage 
supply  requirements,  the  ion  beam  current,  and  thus,  the  detector 
sensitivity  requirements.  Any  constraints  which  are  imposed  on  these 
additional  parameters  strongly  impact  on  the  allowable  values  of  the 
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r-': 


a 

L‘"  , 


i 


yrn 


selectable  parameters  a  ,  R  ,  $  ,  K,  and  d.  Some  guidelines  can  be 

6  6  6 

worked  out  in  advance  without  knowledge  of  any  specific  application. 

The  resolution  goes  to  infinity  as  X2  approaches  zero,  or  as  e2 
approaches  two.  We  are  restricted,  however,  to  values  of  e2  such  that 

X  =  v/| 2  -  e2 1 1  >  |  (5-21) 

Values  of  X  <  £  .,  force  >  2n,  which  is  physically  impossible. 

Also,  as  X  approaches  zero,  the  slits  move  out  toward  infinity  and  we 
again  have  a  size  problem. 

An  examination  of  the  curve  for  l  "la  (see  Figure  20),  shows  the 

z  e 

existence  of  several  poles.  Advantage  can  be  taken  of  the  poles  by 
picking  a  value  of  e2  close  to  the  pole  such  that  l  "  =  .  This 

would  make  the  radial  and  axial  image  planes  coplanar.  Perhaps  the 
best  reason  for  doing  this  is  to  hold  down  axial  extent  of  the  ion 

beam  in  the  radial  image  plane.  This  is  not  essential  but  would  help 

to  avoid  unwanted  losses  at  the  axial  ends  of  the  exit  slit,  which 
must  be  located  in  the  radial  image  plane.  It  also  helps  to  keep  the 
beam  intensity  up  and  to  hold  down  on  the  axial  size  of  the  detector, 
permitting  more  groups  of  detectors  to  be  used. 

A  numerical  search  of  the  region  between  e2  *  2  -  ^  and  the 
first  pole  to  the  left  yields 

e2  »  1.897  (5-22) 

o  "  n  >> 

-Z-  -  —  -  »  3.H6  (5-23) 

ae  ae  ae 

K*R  ,  =9.708  (5-24) 

rel 
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A  numerical  search  of  the  region  between  e  ■  2  +  and  the 
first  pole  to  the  right  yields 

e2  =  2.115  (5-25) 

l  "  £  "  0, 

—  >  —  =  2.951  (5-26) 

ae  ae 

K*R  ,  =  8.708  (5-27) 

rel 

The  first  value  for  e2  [equation  (5-22)]  has  a  higher  resolution 

than  the  second  value  for  e2  [equation  (5-25)].  This  is  at  the 

expense  of  having  the  slits  farther  away,  which  requires  a  larger 

instrument.  It  should  be  remembered  that  these  values  only  pin  down 

e  =  ,  and  not  the  actual  values  of  ae  and  Re. 

Ke 

It  is  the  size  constraints  that  heavily  influence  the  choice  of 
ae  and  Re .  This  breaks  down  into  two  major  considerations:  first, 
how  much  room  is  available  for  the  electrodes;  second,  where  do  the 
slits  have  to  be  with  respect  to  the  electrodes.  If  one  accepts  the 
constraint  imposed  earlier  that  £ "  -  £r"  and  if  the  size  of  the 
electrodes  is  the  primary  concern,  then  one  would  pick  Re,  and 
calculate  ae  from  e2  in  equations  (5-22)  or  (5-25).  If  the  location 
of  the  slits  is  the  primary  concern,  then  one  would  use  equations 
(5-23)  or  (5-26)  to  determine  ag ,  and  then  Rg  would  be  calculated  from 
equations  (5-22)  or  (5-25). 

A  complication  in  choosing  a£  and  Rfi  is  that  the  outer  electrode 

has  an  actual  size  of  r,  ■  a  +  ^  and  R.  ■»  R  +  t  ,  where  d  is  the 

b  e  L  t>  e  i 

electrode  gap.  The  choice  of  d  directly  affects  the  voltage  that 
must  be  applied  to  the  electrodes  to  get  a  given  value  of  Eo,  the 


central  path  electric  field  strength,  via  equation  (2-24).  This 

impacts  on  the  voltage  supply  requirements. 

A  further  complication  is  that  with  a  and  R  chosen  (and  thus  e2), 

e  e 

one  must  now  pick  K  to  get  the  resolution  desired  via  equation  (5-19). 
As  K  gets  smaller,  the  ion  current  goes  down  and  this  impacts  on  the 
detector  sensitivity  required.  The  detector  sensitivity  is  further 
affected  by  the  strength  of  the  ion  source  that  will  be  presented  to 
the  analyzer  entrance  slit. 

The  choices  are  not  easy  and  will  probably  require  many  trials 
before  an  acceptable  compromise  is  reached.  Hopefully,  the  curves  in 
the  Appendix  will  speed  up  the  initial  selection  process. 


F.  Ion  Path  vs  Field  Error 

It  was  initially  assumed  that  all  ions  would  have  trajectories 
near  the  central  path  and  that  the  electric  field  could  be  given  by 
equations  (2-34)  and  (2-35).  However,  as  the  ions  radially  deviate 
from  the  central  path,  the  error  grows  in  the  series  approximation  to 
the  radial  field.  By  comparing  the  series  solution  to  the  exact 
solution  assumed  by  H.  Wollnik  (Ref  4:166)  [equation  (2-20)],  it  is 
possible  to  estimate  the  percent  error  in  the  series  solution  as  a 
function  of  p,  the  deviation  from  the  central  path.  The  percent  error 
is  defined  as 


X  radial  error 


Er(p)  -  foHnik 
EWo¥lnik 


*100 


(5-28) 


with  5  =  0. 


This  equation  has  been  plotted  in  the  Appendix  (see  Figure  26). 
It  should  be  noted  that  the  curve  shows  that  the  percent  error  grows 
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more  rapidly  for  positive  p  than  it  does  for  negative  p.  In  both 
cases,  the  series  solution  underestimates  the  strength  of  the  radial 
field . 

Starting  from  an  on-center  slit,  the  underestimation  of  the  radial 
field  is  going  to  result  in  the  ion  pattern  being  shifted  toward 
lower  values  of  p.  The  pattern  is  fixed  at  p  =  0  because  the  percent 
error  is  zero  for  p  =  0.  Lower  energy  ions  that  traverse  the  p  <  0 
region  will  be  transmitted  in  lower  numbers  than  expected.  Contrawise, 
the  higher  energy  ions  which  traverse  the  p  >  0  region  will  be 
transmitted  in  higher  numbers  than  expected.  In  fact,  due  to  the 
greater  error  when  p  >  0,  more  high  energy  ions  will  be  picked  up  than 
low  energy  ions  are  lost.  This  will  shift  and  distort  the  trans¬ 
mission  curve.  Also,  the  resolution  will  be  slightly  less  sharp  than 
expected  because  of  the  field  errors. 

The  axial  field  series  solution  usually  overestimates  the  true 
field  strength,  but  the  percent  error  is  symmetric  about  C,  *  0,  and 
has  no  effect  on  the  resolution. 

If  one  considers  Figure  12,  one  can  see  that  the  percent  error  in 
the  axial  field  is  given  by 

E  (C)  -  sin@E  (p*) 
z  Wollnik 

Z  axial  error  =  -  .  Q_  -  *100  (5-29) 

sin9E  (p*) 

Wollnik 


where 


* 


P 


H/ 


(1  e2p): 


+ 


(5-30) 


The  overestimation  of  the  axial  means  that  fewer  ions  will  be 


focused  back  onto  the  slit.  As  this  effect  is  weaker  for  p  >  0,  the 
high  energy  ions  will  have  their  transmission  reduced  less  than  the 


lower  energy  ions.  This  increases  Che  distortion  in  the  transmission 
curve  caused  by  the  radial  field  errors.  However,  it  does  nothing  to 
the  resolution. 

A  rough  estimate  of  the  axial  field  errors  can  also  be  gained  by 
the  use  of  equation  (2-8).  This  equation  gives  the  equipotential 
curves  in  the  half  plane  ($  =  constant)  as  curves  of  p  =  f(£).  This 
has  been  plotted  with  e2  set  equal  to  two  (see  Figure  27).  Included 
in  the  graph  is  an  equipotential  curve  through  p  =  0  and  £  =  0,  as 
envisioned  by  H.  Wollnik  (Ref  4:166).  Note  that  the  error  is 
symmetric  about  £  =  0. 

As  the  off-center  slits  are  more  likely  to  direct  their  ions 
through  regions  away  from  the  central  path,  their  resolution  and 
transmission  curves  are  more  likely  to  suffer  distortion  and 
degradation  due  to  the  field  errors  than  are  the  on-center  slits. 

G .  Entrance  Angles  versus  Field  Errors 

If  one  picks  a  maximum  percent  error  that  will  be  tolerated  in  the 

radial  electric  field,  then,  with  the  help  of  Figure  18  in  the 

Appendix,  a  value  for  the  maximum  permissible  p  can  be  chosen.  The 

question  remains:  what  is  the  maximum  permissible  radial  entrance 

I  angle,  «  such  that  p  <  p  for  all  *?  To'  get  an  answer,  we  need 

r  —  max 

the  solution  to  the  radial  equation  of  motion: 

p($)  *  (pi  -yV)  cos (X$)  -  —■  sin(X4>)  +  yj  (3-44) 


By  equation  (4-3),  we  have 


or 


Pi 


(5-31) 


Therefore , 

P(4>)  =  (po - - ^rj  cos(X4>)  -  ~Y~  sin( +  ^y 

'  (5-32) 

The  geometry  of  the  instrument  will  pin  down  £',  a£ ,  and  X.  po 
can  be  set  equal  to  the  radial  center  of  the  entrance  slit  of  interest 
and  n  can  be  calculated  from  equation  (5-1)  for  the  pair  of  slits 
being  considered.  Rho  is  constrained  by  p  <  p  ,  as  selected  from 
Figure  18  in  the  Appendix.  A  computer  search  can  now  be  made,  with  <J> 
as  the  independent  variable,  to  find  the  value  of  '  that  yields 
p  <  Pmax*  Due  to  the  asymmetry  of  Figure  18,  there  exists  two  values 
for  «  ' :  one  for  ions  that  leave  the  entrance  slit  heading  outwards 
(p  >  0),  and  one  for  ions  heading  inwards  (p  <  0). 

By  adding  a  secondary  entrance  slit  behind  the  primary  entrance 
slit,  ions  with  too  large  of  a  radial  entrance  angle  can  be  eliminated 
The  edges  of  the  secondary  slit  are  situated  so  as  to  be  consistent 
with  the  two  values  of  just  computed.  This  will  ensure  that  the 
ions  are  constrained  to  the  region  surrounding  the  central  path  with 
a  percent  field  error  less  than  or  equal  to  the  maximum  error  that  one 
has  chosen  to  tolerate. 

A  similar  calculation  can  be  done  for  the  axial  entry  angle,  «  ' . 
The  solution  to  the  axial  equation  of  motion  is 

CC  _  * 

C,(<p)  *  Cj  cos(  c<f>)  -  — s i n ( g<J) )  (3-51) 
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where 


► 

9  tt 

r 

> 


* 


(4-4) 


or 


(5-33) 


Therefore , 


c(<t>)  = 


■  (“ '  *) 


cos(e<}>) - —  sin(€$) 


(5-34) 


The  geometry  of  the  instrument  will  pin  down  2 ' ,  ag ,  and  e.  j;0 
can  be  set  to  the  axial  center  of  the  entrance  slit,  £0  =  0.  Zeta  is 
constrained  by  L  < L  ,  where  L  is  obtained  from  equation  (5-29). 

A  computer  search  can  now  be  made,  with  $  as  the  independent  variable, 
to  find  the  value  of  '  that  yields  (;  <  Cmax*  Due  t0  the  symmetry 
involved,  only  one  value  of  '  exists. 

The  entrance  slit  is  relatively  narrow  in  the  radial  direction, 
and  the  radial  entrance  angle  permitted  by  the  primary  and  secondary 
slits  does  not  vary  too  much  from  «  '  (max)  for  ions  that  enter  near 
the  radial  edge  of  the  slit.  However,  the  entrance  slit  is  relatively 
long  in  the  axial  direction.  To  hold  down  on  the  deviation  of  the 
maximum  permissible,  axial  entry  angle  from  <*  '  (max),  the  axial 
secondary  slit  edges  need  to  be  far  away  from  the  primary  slit,  i.e., 
as  close  to  the  electrodes  as  possible.  This  can  be  accomplished  with 
a  single  secondary  slit  located  near  the  electrodes  with  both  the 
appropriate  axial  and  radial  edges. 

The  existence  of  the  secondary  slit  accomplishes  four  things: 
first,  it  limits  ion  trajectories  to  the  central  region;  second,  it 
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prevents  ions  from  straying  from  one  group  of  exit  slits  to  another 
group  causing  an  ambiguity  in  the  data;  third,  it  defines  the  soli 
angle  of  acceptance  of  the  primary  entrance  slit;  and  fourth,  it 
serves  to  hold  down  on  the  fringe  fields.  It  might  also  be  desirable 
to  place  a  matching  secondary  slit  on  the  exit  end  of  the  electrodes. 

H .  Total  Beam  Current 

The  total  beam  current  must  be  arrived  at  with  statistical 
mechanics.  The  procedure  will  be  qualitatively  outlined  below.  No 
exact  equations  will  be  given  since  the  problem  is  too  dependent  on 
the  particular  experiment  for  any  one  set  of  equations  to  be  of  much 
use.  There  are  four  statistical  populations  that  must  be  worked  out 
in  order  to  get  the  beam  current. 

The  first  population  is  the  ionized  plasma  that  is  under 
observation.  This  part  of  the  problem  is  solely  a  function  of  the 
experiment  being  conducted. 

The  second  population  is  the  plasma  that  comes  through  the  primary 
slit.  This  part  of  the  problem  is  essentially  that  of  a  gas  escaping 
through  a  small  hole  in  a  container  with  an  unlimited  supply  of  gas  in 
the  container.  One  must  also  take  into  account  whether  or  not  there 
is  a  net  motion  of  the  plasma  relative  to  the  primary  entrance  slit. 

If  the  plasma  is  dense  or  hot,  as  in  a  pulsed  plasma  device,  then 
plasma  effects  such  as  magnetohydrodynamics  may  have  to  be  considered. 
In  general,  the  second  population  will  have  a  higher  average  energy 
than  the  first  population.  It  will  also  have  an  angular  distribution 
that  is  peaked  in  the  forward  direction. 

The  third  population  is  that  portion  of  the  second  population 
that  clears  the  secondary  slit.  It  can  be  determined  by  integrating 


the  second  population  over  the  solid  angle  defined  by  the  primary 
and  secondary  slits. 

The  fourth  population  is  that  portion  of  the  third  population 
that  clears  the  exit  slit.  It  can  be  determined  by  multiplying  the 
third  population's  distribution  function  by  the  transmission  function 
and  then  integrating  over  energy.  The  result  of  the  integration  is 
the  total  beam  current  as  seen  by  the  detector. 

If  the  primary  slit  width  and  length  have  been  carried  along  as 
unspecified  parameters,  then  we  will  now  have  the  beam  current  as  a 
function  of  the  slit's  width  and  length.  It  will  now  be  possible  to 
make  the  trade-off  between  detector  sensitivity  and  the  resolution  by 
selecting  the  appropriate  slit  width  and  length. 

The  only  problem  with  all  of  this  is  that,  if  the  first 
population's  distribution  function  were  known,  then  we  would  not  need 

rr 

the  analy  •  .v  at  ail.  If  an  approximate  model  is  available,  the  above 
procedure  can  be  used  to  select  the  primary  slit.  In  general,  one 
needs  to  start  with  the  detector  and  work  backwards  so  as  to  connect 
the  detector  data  to  the  original  population.  This  can  be  made 
simpler  if  one  assumes  that  the  original  population  obeys  Maxwell- 
Boltzmann  statistics.  It  is  still  essential  to  know  how  the  plasma 
is  moving  with  respect  to  the  primary  slit. 

One  fact  that  simplifies  things  is  that  the  resolution  and 
transmission  are  in  terms  of  the  ion  kinetic  energy,  T,  outside  the 
field  and  the  kinetic  energy,  T0 ,  of  an  ion  that  obeys  the  centripetal 
force  equation.  This  gets  one  outside  of  the  electrodes  but  inside 
the  secondary  entrance  slit.  The  secondary  entrance  slit  only  clips 
the  angular  distribution  and  does  not  really  affect  the  energy 
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spectra.  Thus,  the  detector  data  can  be  directly  related  to  the 
number  of  ions  passing  through  the  primary  slit  with  a  specific  energy. 

I .  Conclusion 

A  toroidal  geometry,  electrostatic  analyzer  has  several  distinct 
advantages:  first,  it  can  be  used  as  a  multi-channel  analyzer 

sampling  different  segments  of  the  energy  spectra  at  the  same  time; 
second,  the  range  of  energies  that  can  be  sampled  simultaneously  is 
more  than  an  order  of  magnitude  wider  than  it  is  with  either  a 
cylindrical  or  a  spherical  system;  and  third,  the  resolution  is 
significantly  higher  than  it  is  with  either  a  cylindrical  or  a 
spherical  system.  These  advantages  definitely  enhance  the  speed  and 
accuracy  of  the  electrostatic  analyzer. 

J .  Recommendations 

My  recommendations  are  twofold.  First,  the  fringe  fields  need  to 
be  treated  more  accurately.  Perhaps  Albrecht's  work  (Ref:3) 
will  be  of  some  help.  Second,  the  image  aberrations  need  further 
studying.  This  will  necessitate  an  investigation  of  higher  order 
terms  in  the  electrostatic  equations,  the  mechanic  equations,  and  the 
optic  equations. 

/ 


74 


1.  Svartholm,  N.  "Velocity  and  Two-Directional  Focusing  of  Charged 
Particles  in  Crossed  Electric  and  Magnetic  Fields,  I," 

ARKIV  FOR  FYSIK,  Band  2  nr  20;195-207  (1950). 

2.  Ewald,  H.  and  Liebl,  H.  "Der  Astigraatismus  der  Toroid-Kondensators 
2.  Naturforsch,  10a:872-876  (1955). 

3.  Albrecht,  R.  "Das  Potential  in  doppelt  geKrummten  Kondensatoren ," 

2  .  Naturforsch,  lla:156-163  (1956). 


Wollnik,  H.  "Electrostatic  Prisms,"  Focusing  of  Charged  Particles, 
Volume  2,  edited  by  A.  Septier.  New  York:  Academic  Press,  1967. 

Duckworth,  H.  E.  and  Ghoshal,  S.  N.  "High-Resolution  Mass 
Spectroscopes,"  Mass  Spectrometry,  edited  by  C.  A.  McDowell. 

New  York:  McGraw  Hill  Book  Company,  Inc.,  1967. 

Dahl,  P.  Introduction  to  Electron  and  Ion  Optics,  New  York: 
Academic  Press,  1973. 

Paolini,  F.  R.  and  Theodoridis,  G.  C.  "Charged  Particle  Trans 
mission  through  Spherical  Plate  Electrostatic  Analyzers," 

The  Review  of  Scientific  Instruments,  38:579-588  (May,  1967). 


Shown  on  the  following  pages  are  graphs  which  should  help  the 

reader  to  design  his  own  instrument.  Most  of  these  curves  have  been 

plotted  with  a  common  e2-axis,  where  e2  is  the  ratio  of  the  two  radii, 

,  of  the  mid-electrode  surface.  It  should  be  noted  that  several  of 
Ke 

the  curves  have  been  divided  by  ae ,  to  get  curves  that  are  of  greater 
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Figure  13  -  Upper  Energy  Limit  of  Off-Center  Slits 


f  DEL.  TR 


T  \/(  E*UG 


Figure  15  -  Energy  Range  of  Off-Center  Slits 


Figure  17  -  Location  of  Axial  Focal  Point 
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Figure  19  -  Location  of  Radial  Image  Plane 
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Figure  20  -  Location  of  Axial  Image  Plane 
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Figure  21  -  Radial  "Focal  Length 
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Figure  24  -  Resolution 
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